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■Q^"; ABSTRACT 

^^ • We present BeppoSAX follow-up observations of GRB980425 obtained with the Narrow Field Instru- 

ments (NFI) in April, May, and November 1998. The first NFI observation has detected within the 8' 
radius error box of the GRB an X-ray source positionally consistent with the supernova 1998bw, which 
exploded within a day of GRB980425, and a fainter X-ray source, not consistent with the position of the 
supernova. The former source is detected in the following NFI pointings and exhibits a decline of a factor 
of two in six months. If it is associated with SN 1998bw, this is the first detection of X-ray emission 

CO ' from a Type I supernova above 2 keV. The latter source exhibits only marginally significant variability. 

The X-ray spectra and variability of the supernova are compared with thermal and non-thermal models 
of supernova high energy emission. Based on the BeppoSAX data, it is not possible to firmly establish 
, which of the two detected sources is the GRB X-ray counterpart, although probability considerations 

*f^ , favor the supernova. 

rf-\ ' Subject headings: gamma rays: bursts — Supernovae: individual: SN 1998bw 
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O I 1. INTRODUCTION with the "classical" firebaU model (e.g., Rees & Meszaros 

5^ ■ The GRB of 1998 April 25, detected both by the Bep- 1992, Piran 1999) GRB980425 has been regarded as a 
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Tj'oi^ lo/f/^oo r n nnoc rp- , , i nno\ rpi, bascd On individual cascs of supcrnovae With Outstanding 

ESO 184-G82 [z = 0.0085, Imney et ai. 1998). Ihe ,. , ; • /^ 7 , 1r^r^r^ rr^ ;; ; r 

, 1 , .,.; r u • -J u J. /i i optical properties (Germany et al. 1999; luratto et al. 

low probability of a chance coincidence between the two ^^_^. r^il, .^. , ^„}. „ „ ri.,_,,„ -,nnn^ 



, / iri-4 r< ^ i i inno\ T, i -ti, j 4-1, 1999; Terlevich, Fabian, & Turatto 1999). 
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rrt ' hypothesis 01 a physical association between GRB980425 ., -i^rr^^ ° , •;! ,, t-> r. a ^^ at tt- , , t 

■ ■ ' d th WFC error box with the BeppoSAX Narrow Field In- 

TT ■ ii, iu riTKTD c u- u J ^,■fi- struments (NFI) were immediately activated, starting 10 

However, since the other GRBs tor which a redshitt , . ^ , n . i . mi i ; ; • r . 

, . -1 i„i 1 i J i 1 J- i hours and one week after the event. Ihe detection ot two 

measurement is available are located at larger distances -11 ^r r 1 • 1 i • 

/ r,'7 ^■u^ J 1, i-ju 1 J previously unknown X-ray sources - one ot which being 
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GRB X-ray counterpart, and, as a secondary though not 
less important scope, at monitoring the X-ray emission of 
SN 1998bw, considering the pecuharity of this object and 
the poor knowledge of the X-ray behavior in supernovae 
in general (see review by Schlegel 1995), and particularly 
of Type I. 

In this Letter we present the high energy characteris- 
tics of the prompt event as measured by the BeppoSAX 
GRBM and WFC, and the results of the follow-up NFI 
observations (§2) and discuss their implications in view of 
the detection of SN 1998bw in the GRB field (§3). A pre- 
liminary report on these data has been given in Pian et al. 
(1999). 

2. DATA ANALYSIS AND RESULTS 
2.1. Prompt event 

GRB980425 triggered the BeppoSAX GRBM at 
21:49:11 UT, and was simultaneously detected by the Bep- 
poSAX WFC unit 2 (Sofhtta ct al. 1998). The event had 
a duration of 31 s in the range 40-700 keV and 40 s in the 
range 2-26 keV. It exhibited a single, non structured peak 
profile in both bands, with the indication of a ~5 s soft lag 
(Fig. 1). Some flux brightening and successive decrease, 
lasting altogether ~^10 seconds, are seen in the WFC light 
curve after the first 40 seconds, but not in the 7-rays. 

The GRBM and WFC light curves appear well corre- 
lated, with the indication of a ~5 seconds lag of the lower 
energies with respect to the higher energies. Figure 2 
shows the Discrete Correlation Function (DCF) between 
the two light curves. This correlation method is suited to 
search for correlations and temporal lags between two dis- 
crete, and possibly unevenly sampled, data trains (Edel- 
son & Krolik 1988). To a maximum of the DCF amplitude 
at a positive temporal lag corresponds a positive correla- 
tion, with the higher energies leading the lower ones. A 
maximum at a positive temporal lag of ~5 seconds is evi- 
dent. The asymmetric shape of the DCF function ampli- 
tude around its peak reflects the fact that the flux decay 
after maximum is slower at X-rays than at 7-rays (see Fig. 

!)■ 

The fluences are (2.8 ± 0.5) xlG ^ erg cm ^ and 
(1.8 ± 0.3) xlO-s erg cm'^ in the 40-700 keV and 2-26 
keV energy range, respectively. (The Galactic absorp- 
tion in the direction of GRB980425, Nhi ^ 4 x 10^° 
cm^-^, from Schlegel, Finkbeiner, & Davis 1998, is neg- 
ligible at energies higher than 2 keV.) The spectral index 
of the 7-ray spectrum, averaged over the burst duration, is 
a = 1.2 ± 0.2 (/i, (X i^~"), and that of the X-ray spectrum 
is a = 0.41±0.25. Strong spectral softening is evident dur- 
ing the event (Frontera et al. 1999). No peculiar temporal 
or spectral characteristics are noted in this GRB. 

The spacecraft aspect reconstruction conditions allowed 
us to only poorly constrain the WFC error box of the GRB, 
which has a radius of 8'. This made the search of an X- 
ray transient more difficult. The IPN allowed a substantial 
reduction of the GRB error box (see Galama et al. 1999b). 

2.2. NFI Target of Opportunity Observations 

The BeppoSAX NFI were pointed at the 8' radius error 
box determined by the WFC at three epochs starting 10 
hours after the GRB: April 26-28, May 2-3 and November 
10-12 (see Journal of Observations in Table 1; note that 



the April pointing was uninterrupted, but has been split in 
two parts only for the data analysis purpose). The strategy 
of performing NFI pointings at so large time intervals af- 
ter the GRB (one week and six months), besides promptly 
thereafter, is not usually adopted for GRBs and was dic- 
tated by the ambiguity of the results obtained with the 
April observation, which suggested a quite different case 
than previously observed X-ray afterglows. 

Event files for the LEGS and MEGS experiments were 
linearized and cleaned with SAXDAS at the BeppoSAX 
Science Data Center (SDC; Giommi & Fiore 1998). LEGS 
(0.1-4 keV) and MEGS (unit 2 and 3, 1.6-10 keV) images 
for each pointing were extracted using the XIMAGE soft- 
ware package. The MEGS images, better exposed and of 
higher signal-to-noise ratio than the LEGS ones, are re- 
ported in Fig. 3. 

The analysis of the MEGS imaging data of the first por- 
tion of the first pointing (Fig. 3a) shows that inside the 
intersection of the WFC error circle and IPN annulus, two 
point-like, previously unknown X-ray sources are detected 
with a positional uncertainty of 1'.5: 1SAXJ1935. 0-5248 
(hereafter SI), at RA ^ 19h 35m 05.9s and Dec = -52° 
50' 03", and 1SAXJ1935.3-5252 (hereafter S2), at RA = 
19h 35m 22.9s and Dec = -52° 53' 49". Note that the 
coordinates distributed by Pian et al. (1998) have been 
revised in November 1998, to take into account a system- 
atic error due to the non-optimal spacecraft attitude dur- 
ing the April and May 1998 observations (see Piro et al. 
1998a). The revised position of SI is consistent within the 
uncertainty with the position of SN 1998bw detected in 
the WFC error box (Galama et al. 1998; Kulkarni et al. 
1998a) and exploded simultaneously with the GRB with 
an uncertainty of ~1 day (Iwamoto et al. 1998), while the 
revised position of S2 is ~ 4'. 5 away from SN 1998bw, and 
therefore inconsistent with it (see Fig. 1 in Galama et al. 
1999b). 

The LEGS and MEGS count rates and upper limits for 
both sources during the three pointings have been com- 
puted within circles of 3' radius and corrected for the 
local background estimated within circles of similar size 
(Table 1). The upper limits have been estimated as ex- 
plained in the Appendix. This method takes into account, 
in addition to the normal photon statistics, also the fact 
that, at these flux levels, the LEGS and MEGS background 
may be dominated by the fluctuations of the cosmic X-ray 
background. 

Source SI is detected by the MEGS also in the follow- 
ing pointings at a position consistent with that of the first 
observation. The observation of November 1998 shows a 
decrease in the X-ray flux of approximately a factor of two 
with respect to the level measured in April-May 1998. The 
source is also detected by the LEGS in the April pointing, 
and not in the November pointing (LEGS data of the May 
1998 pointing are excluded from the analysis because of 
the extremely short exposure time). 

Due to the limited spatial resolution of the LEGS and 
MEGS detectors and to the faint emission level of source 
S2, estimating its flux is made difficult by the background 
contamination and by the proximity of the brighter source 
SI. During the second portion of the April pointing, as 
well as in the November 1998 pointing, S2 is not detected 
by the MEGS, while signal from a position consistent with 
that of S2 in April is marginally detected in the May 1998 



pointing (Fig. 3c), with lower flux than in the first obser- 
vation (see Table 1), but consistent with it at the ^2-0- 
level. We note that the May detection, albeit of a signal- 
to-noise formally lower than 3, has a very low probability 
of being a background fluctuation (~ 10~^) when consid- 
ered together with the April detection. However, we con- 
servatively report in Table 1 also the 3-a upper limit for 
the May measurement. All upper limits to the flux of S2 
are consistent with the level of the detection. There is no 
significant detection of S2 in the individual LEGS images. 

No significant signal above background is detected 
by the BeppoSAX high energy instruments PDS and 
HPGSPC in any of the three pointings. 

Light curves and spectra for each pointing were accumu- 
lated with the XSELECT tool, using a 3 arcmin extraction 
radius both for the LEGS and the MEGS, which provides 
only 80% of flux, but allows partially avoiding the mutual 
contamination of SI and S2. Since the local background 
intensity is similar to that measured from files accumu- 
lated from blank fields available at the SDG, we used the 
latter, which are affected by a smaller uncertainty. No 
significant variability within each pointing is exhibited by 
either source. 

Spectral analysis of the LEGS and MEGS data has been 
done with the XSPEG 10.0 package using the response 
matrices and auxiliary files available at the SDG. LEGS 
and MEGS spectra flux distributions of both sources have 
been grouped in order to achieve a signal-to-noise ratio of 
at least ~3 in each bin. 

For SI, a fit with a single power-law (F^ ex ly^") ab- 
sorbed by Galactic extinction of the individual MEGS 
spectra is satisfactory, with reduced x^ values well below 
one, due to the large errors. For these fits at individual 
epochs the LEGS spectra have not been used because they 
have too low signal-to-noise ratio. 

Since no significant spectral variability is seen from 
epoch to epoch, we averaged the LEGS and MEGS spectra 
of April and May, to increase the signal-to-noise ratio, and 
fitted a single power-law to the average rebinned spectrum. 
We excluded the November spectrum from this average be- 
cause the count rates indicate that the flux varied at that 
epoch with respect to April-May (Table 1). 

We obtain a spectral index of a = 1.0 ± 0.3 (all fit pa- 
rameters for the average spectrum are reported in Table 2). 
The fitted Nhi is found to be consistent with the Galactic 
value and therefore we have fixed it to that value. The 
intercalibration constant between the LEGS and MEGS 
data is within the expected range. However, the fit is for- 
mally not completely satisfactory (reduced x^ = 1-2, see 
Table 2), due to a flux excess at energies below ~ 1 keV 
(Fig. 4a). ^ ^ 

We also tried a fit with a thermal bremsstrahlung with 
Galactic absorption, obtaining a temperature kT ~ 8 keV 
and a not satisfactory reduced x^ of 1.7, still due to the 
presence of a soft excess (Fig 5a). 

Fitting the data assuming no Galactic absorption still 
yields some soft excess with respect to both a power-law 
and a bremsstrahlung model, therefore we tend to exclude 
that the effect is due to an overestimate of the Galactic 
extinction in the direction of SN 1998bw (as possibly sug- 
gested by Patat et al. 1999). 

A fit of the data with a broken power-law and Galactic 
absorption yields an index ai = 2.0 ± 0.5 below 1.4 keV 



and a2 — 0.7 ± 0.4 at the higher energies. The reduced 
X^ is 0.5, considerably lower than for the single power-law 
and for the bremsstrahlung models (see Fig. 6). 

To account for the soft excess we also tried composite 
fits of a power-law or a thermal bremsstrahlung with a 
black body model plus Galactic absorption. The former 
fit yields a spectral index a = 0.5 ± 0.3 and a black body 
temperature kT = 90 ± 20 eV, with a reduced x^ = 0.5 
(Fig. 4b). The latter yields a bremsstrahlung temperature 
of kT ~ 16 keV (see Table 2), and black body tempera- 
ture and normalization consistent with those found in the 
former case (reduced x^ = 0.5, Fig. 5b). 

We found no obvious reason for the excess to be spuri- 
ous, such as instrumental effect or non optimal background 
subtraction. Therefore we considered it real, although not 
highly significant (see Figures 4a and 5a). The chance 
probabilities that adding a black body component either 
to a power-law or to a thermal bremsstrahlung model im- 
proves the fit are 1.3% and 0.17%, respectively. Fitting 
the data with a broken instead of a single power-law has 
a chance probability of improvement of 4%. 

The low signal-to-noise ratio of the LEGS data in 
November 1998 does not allow us to test the goodness of 
the fit of those data with a black body model. By fixing the 
black body temperature to the best fit value of the April- 
May spectrum, 90 eV, and fitting a power-law plus black 
body model to the November MEGS spectrum and LEGS 
upper limit we get a black body normalization upper limit 
consistent with the value obtained for April-May, indicat- 
ing that this component has not varied significantly. The 
power-law normalization has instead varied by a factor of 
-2. 

Although there is no detection of S2 in the individual 
LEGS images, signal is present in the April-May coadded 
image. Therefore, since the MEGS spectra of the individ- 
ual observations have a very low signal-to-noise ratio, we 
only used the rebinned LEGS-I-MEGS spectrum obtained 
from the average April-May image, and fitted it with a 
power-law of index a = 1.5 ± 0.4 plus Galactic absorption 
(x' = 0.5). 

Due to the high uncertainty in the power-law spectral 
index and to the low level of both SI and S2, we estimated 
their intensities in the 2-10 keV band by assuming a stan- 
dard factor of conversion from count rates of 9.3 x 10"^^, 
corresponding to an adopted power-law spectral shape of 
index 0.5. The l-a uncertainties on the intensities have 
been obtained by similarly scaling the errors on the count 
rates (Table 1). These intensities are reported in Figure 7, 
along with the WFG light curve in the 2-10 keV band, ob- 
tained by binning in intervals of 5-20 seconds the temporal 
profile reported in Figure 1. Differences between these in- 
tensities and those obtained by adopting a bremsstrahlung 
model do not exceed 10%. Note that the contribution of 
the black body component to the intensity of SI in the 
2-10 keV range is negligible. 

3. DISCUSSION 

3.1. Light curve and spectral shape of source SI 

The X-ray light curve of source SI measured by the NFIs 
shows a decay of a factor of two in ~6 months (Fig. 7b) , 
much slower than X-ray GRB afterglows so far observed. 
Assuming, as suggested by the positional coincidence and 



by variability, that SI is associated with SN 1998bw, this is 
the first detection of medium energy X-ray emission from 
a Type I supernova (there is a unique case of Type I su- 
pernova detected in soft X-rays, the Type Ic SN 19941, 
Immler et al. 1998a) and the earliest detection of X-rays 
after supernova explosion. 

At the distance of SN 1998bw, 38 Mpc, the luminos- 
ity observed in the range 2-10 keV, ~ 4 — 7 x 10'*" erg 
s~^, would be compatible with that of other supernovae 
detected in the same energy band (Kohmura et al. 1994; 
Houck et al. 1998; Schlegel 1995, and references therein). 

However, the observed luminosity and variation thereof 
represent only an upper limit to the luminosity and a lower 
limit to the amplitude of X-ray variability of SN 1998bw, 
respectively, due to the possible contribution of its host 
galaxy, a face-on spiral galaxy about one tenth of the size 
of our Galaxy, which is only very marginally resolved in 
the BeppoSAX data. In fact, a galaxy of that type and 
size could easily account for almost all of the X-ray emis- 
sion observed in November 1998, when the flux was lowest 
(see e.g., Fabbiano 1989). 

The observed decay of SI in the 2-10 keV band is well 
fitted by a power-law F{t) oc t'P with p = 0.16 ± 0.04 
(reduced x^ — 0.7). The fit with an exponential law 
F{t) ex e~*/'^ with (3 = 500 ± 100 days has a reduced x^ of 
2.4, corresponding to a probability of ^10% for 2 degrees 
of freedom, therefore not negligible (Fig. 7b). Both trends 
would be similar, considering the unknown dilution by the 
host galaxy of SN 1998bw, to the X-ray behavior of other 
supernovae (e.g., Kohmura et al. 1994; Zimmermann et 
al. 1994; Houck et al. 1998) and predicted by models of 
thermal bremsstrahlung of energetic electrons within the 
circumstellar medium (see e.g., Chevalier & Fransson 1994; 
Chugai & Danziger 1994). 

The prompt X-ray emission observed for SN 1998bw 
requires that the circumstellar medium is highly ionized 
(perhaps by the powerful explosion), to allow the X-rays to 
escape so soon after the explosion (see Zimmermann et al. 
1994), and also very dense, as inferred also from the large 
radio output (Kulkarni et al. 1998a; Wieringa, Kulka- 
rni, & Frail 1999). In these conditions, the X-rays might 
be produced by the reverse shock which results from the 
pressure of swept-up material on the outgoing shock and 
propagates back into the shocked supernova gas (Chevalier 
& Fransson 1994; Schlegel 1995). 

The temperature obtained from the thermal 
bremsstrahlung fit to the BeppoSAX LECS-hMECS spec- 
tra, admittedly very poorly constrained, is compatible 
with that fitted to the medium or hard X-ray spectra of 
other supernovae (Kohmura et al. 1994; Leising et al. 
1994; Dotani et al. 1987). 

On the other hand, the mildly relativistic conditions 
evidently present in the expanding shock of SN 1998bw 
at early epochs (Kulkarni et al. 1998a) and the accept- 
able power-law fit obtained for the medium energy X-ray 
spectra might suggest that non-thermal mechanisms are 
responsible for the X-ray emission such as synchrotron ra- 
diation by the extremely energetic electrons, as was mod- 
eled for the radio emission by Li and Chevalier (1999), 
or inverse Compton scattering of relativistic electrons off 
optical/UV photons of the thermal ejecta (see Canizares 
et al. 1982). The X-ray spectral index is consistent with 



that measured for the radio spectrum starting ^15 days 
after the explosion (Kulkarni et al. 1998a; Wieringa et al. 
1999), and with the slope connecting quasi-simultaneous 
radio and X-ray measurements (a ~ 0.8). Therefore, it 
is difficult to establish whether the X-rays are produced 
through synchrotron or inverse Compton radiation. How- 
ever, if inverse Compton losses were dominant, radio emis- 
sion production would be rapidly inhibited (Schlegel 1995), 
contrary to the observations. 

Assuming the X-rays have a synchrotron origin and 
adopting a single power-law of index a ~ 0.8 for the radio- 
to-X-ray synchrotron spectrum, we obtain a bolometric lu- 
minosity of ^ 9 X lO**" erg s"-'^. Assuming that the radio 
and X-ray emitting regions are cospatial and expanding 
with a speed of ~ 0.3c (Kulkarni et al. 1998a), we de- 
rive a magnetic field of ^ 1 Gauss, similar to that found 
for SN 1980K by Canizares et al. (1982) assuming inverse 
Compton losses are responsible for the X-ray production. 

However, the limited signal-to-noise ratio of the Bep- 
poSAX spectra does not allow us to choose between syn- 
chrotron radiation and thermal bremsstrahlung as the 
mechanism for the medium energy X-ray production. 

The emission component detected in the softer part of 
the BeppoSAX spectrum, if fitted with a black body, has 
a temperature of ~0.1 keV, corresponding to a black body 
total luminosity of ^ 10*** erg s~^. The inferred linear 
size of the emitting region is about one third of the solar 
radius, too large for a compact object left as a remnant 
of the supernova explosion, but approximately compati- 
ble with the size of the putative accretion disk promptly 
formed as a consequence of the "hypernova" or "coUapsar" 
phenomenon, of which SN 1998bw might be an example 
(Paczyhski 1998; MacFadyen & Woosley 1998; Woosley, 
Eastman, & Schmidt 1999a). However, such a compact 
object or disk could be hardly visible at a so early epoch, 
due to the optical thickness of the material enshrouding 
it. 

Black body soft X-ray emission is not expected from 
the supernova itself or from the expanding shell or ejecta, 
due to non-equilibrium conditions of the system. How- 
ever, the occurrence of short (^1000 s) thermal bursts at 
early times after supernova explosion has been predicted, 
albeit at very soft energies, possibly even lower than those 
observable by the LEGS (Schlegel 1995). We found no ev- 
idence in the LEGS light curves of similar events, but the 
sampling is not conducive to that detection. 

However, since there is no evidence of variability of the 
soft component in the long term, it might not be related to 
the supernova. In fact, its spectrum can be fitted also with 
a power-law superimposed on, and steeper than, the one 
which describes the spectrum at higher energies. This sug- 
gests that the component might have a more complex spec- 
trum (possibly extending toward ultraviolet wavelengths), 
of which a black body or steep power-law are only approxi- 
mations. It could rather be a persistent (or slowly variable 
with a small amplitude) source of soft X-rays such as the 
host galaxy itself, or just its bulge, or the superposition 
of unresolved X-ray sources within that galaxy, or diffuse 
hot gas, or the HII region in which SN 1998bw is located 
(Galama et al. 1998), or the underlying cluster DN 1931- 
529, or more probably, the sum of some or all of these con- 
tributions. The limited angular resolution of BeppoSAX 
does not allow us to disentangle this component from the 



supernova itself. Notwithstanding this possibihty, we note 
that the unabsorbed luminosity in the 0.1-2 keV range, 
5 X 10"*° erg s""'^, given by the superposition of the fitted 
power-law and black body components, is similar to lumi- 
nosities of supernovae observed in soft X-rays (Canizares 
et al. 1982; Bregman & Pildis 1992; Zimmermann et al. 
1994; Schlegel et al. 1996; Fabian & Terlevich 1996; Imm- 
ler et al. 1998b). 

3.2. The association between GRB980425 and 
SN 1998bw 

The GRB980425 prompt event is relatively weak with 
respect to other GRBs, and rather soft. However, it has 
no outstanding features with respect to other BeppoSAX 
or BATSE GRBs, which might suggest a peculiar counter- 
part at longer wavelengths, such as a bright supernova, in- 
stead of a "classical" power-law fading afterglow. The ^5- 
seconds temporal lag between the WFC and the GRBM 
light curves could be due to a delay of X-ray emission 
during the burst with respect to the 7-rays, or ascribed 
to intrinsic absorption in a medium becoming increasingly 
transparent (see e.g., Bottcher et al. 1999). Similar soft 
lags from few to ^10 seconds are observed also in other 
GRBs (Piro et al. 1998b; Piro et al. 1998c; Frontera et al. 
1999). 

If SN 1998bw is the counterpart of GRB980425, the pro- 
duction of 7-rays could be accounted for by the explosion 
of the 14 Mq helium core of a ~ 35 Mq star (Woosley et 
al. 1999a; MacFadyen & Woosley 1999) and by the sub- 
sequent expansion of a relativistic shock, in which non- 
thermal electrons are radiating photons of ^100 keV, pro- 
vided the explosion is asymmetric, i.e. the GRB is pro- 
duced in a relativistic jet (Iwamoto et al. 1998; Woosley 
et al. 1999a; Hoflich, Wheeler, & Wang 1999; Rej 1998; see 
however, Kulkarni et al. 1998a). The presence of an unde- 
tectable, or barely detectable, non-thermal GRB remnant, 
underlying the brighter thermal supernova ejecta cannot 
be excluded (see e.g., Iwamoto 1999). 

Recent speculations have led to the proposal that every 
long (> 1 s) GRB is formed via supernova, or hypernova, 
explosion (MacFadyen & Woosley 1999). The presence of 
a supernova underlying the GRB afterglow has been re- 
cently tested for the optical transients of some GRBs with 
suggestive results (GRB970228, Reichart 1999; Galama et 
al. 1999c; GRB970508, Germany et al. 1999; GRB980326, 
Bloom et al. 1999; GRB990510, Fruchter et al. 1999d; 
Beuermann et al. 1999; GRB990712, Hjorth et al. 1999). 
Indeed, the recent discovery of a GRB optical counter- 
part at the intermediate redshift z = 0.43 (Galama et 
al. 1999d) might support a continuity of properties be- 
tween GRB980425 and the other precisely locahzed GRBs, 
perhaps based on the different amount of jet coUimation 
(Woosley, MacFadyen, & Heger 1999b) or different beam- 
ing, depending on the degree of jet alignment (Eichler & 
Levinson 1999; Cen 1998; Postnov, Prokhorov, & Lipunov 
1999). In highly coUimated or highly beamed GRBs the 
non-thermal multiwavelength afterglow could overwhelm 
the underlying supernova emission. This should instead 
be detected more clearly in less collimated or less beamed 
(i.e., seen off-axis) GRBs, like GRB980425, which are, or 
appear, weaker. Assuming association with SN 1998bw 
and isotropic emission, the total energy of GRB980425 in 
the 40-700 keV, ~ 5 x 10^^ erg, is at least four orders of 



magnitude less than that of GRBs with known distance. 

On the other hand, disregarding the fact that the 
probability of a chance coincidence of GRB980425 and 
SN 1998bw is extremely low, one might consider S2 as 
the X-ray counterpart candidate of the burst. The possi- 
ble detection of S2 in 2-3 May 1998, one week after the 
GRB, implies, with respect to the first NFI detection in 
April, a much slower decay than that normally observed 
for X-ray afterglows (e.g., Costa et al. 1997; Nicastro et 
al. 1998; Dal Fiume et al. 1999; in 't Zand et al. 1998; 
Nicastro et al. 1999; Vreeswijk et al. 1999; Heise et al. 
1999). 

Assuming a power-law decay between the X-ray flux 
measured by the WFC in the 2-10 keV range in the last 
~20 seconds of the GRB and the flux measured in the 
flrst NFI observation (Fig. 7a), we derive a power-law 
index p ~ 1.5, which is similar to commonly observed X- 
ray afterglows. If S2 is an afterglow, one would expect 
that its intraday variability followed this same temporal 
behavior. Therefore, we have binned the light curve of S2 
in the first portion of the April 1998 pointing in five inter- 
vals of 20000 seconds each. We have then connected with 
power-laws the last WFC measurement with the first and 
last of these fluxes and have determined their indices to 
be p ~ 1.6 and p ~ 1.4 (Fig. 7c), respectively. The re- 
duced x^ values computed for these two power-laws with 
respect to the remainder four NFI data points of April are 
3 and 30, respectively, corresponding to low probabilities 
(1% and <Cl) that the power-laws describe the observed 
intrapointing light curve. (All points seem rather consis- 
tent with a constant trend.) 

The upper limit derived for the second portion of the 
April pointing is inconclusive. However, the detection of 
S2 in May 1998 suggests a marginal deviation from the 
above power-laws ( ^ 2.5a). Therefore, the present data 
exclude at a confldence level of ^99%, or higher, that S2 
is an afterglow, unless a small re-bursting, one week after 
the GRB, is superimposed to the power-law monotonic de- 
cline. This would be reminiscent of GRB970508, although 
the time scales for re-bursting occurrence and duration 
would be very different (Piro et al. 1998b). 

4. CONCLUSION 

Two previously unknown X-ray sources have been de- 
tected by the BeppoSAX NFIs in the fleld of GRB980425. 
Neither of them has the obvious characteristics of an X-ray 
afterglow, when compared to previously observed cases. 

SN 1998bw is a very interesting candidate for further 
monitoring in the X-rays. Thanks to its rapid slew capabil- 
ity and to its wide energy range, BeppoSAX has promptly 
measured its spectrum up to energies beyond ^ 5 keV, 
where supernovae have been so far largely unexplored. Un- 
fortunately, the signal is relatively modest (this is the sec- 
ond most distant supernova detected so far in the X-rays, 
after SN 1988Z, located at 95 Mpc, Fabian & Terlevich 
1996), and therefore longer exposures and data of better 
signal-to-noise ratio are necessary to study in detail this 
source and its environment. 

Concerning the identification of the X-ray counterpart of 
GRB980425, our tentative conclusion is that SI has a high 
probability of being associated with GRB980425, while S2 
is more probably a variable fleld source, albeit constant in 
the long term, like an active galaxy or a Galactic X-ray bi- 



nary (we note that the probability of detecting by chance 
a source of the level of S2 is rather high, ^10%, Cagnoni 
et al. 1998; Giommi et al. 1998). A spectroscopic survey 
of the NFI error box of S2 has been inconclusive in this 
respect (Halpern 1998). Observations of this field by an 
X-ray instrument with higher sensitivity (e.g., XMM) and 
spatial resolution (e.g., Chandra) than those attained by 
BeppoSAX might help elucidating this controversial issue. 

We thank the BeppoSAX Mission Planning Team and 
the BeppoSAX SDC and SOC personnel for help and sup- 
port in the accomplishment of this project. 



APPENDIX 

ESTIMATE OF UPPER LIMITS AT 3-cr CONFIDENCE LEVEL 

Suppose that m counts are found in the region where the target is expected to appear and that in the same search area 
b counts are expected from the background. If m — 6 < 3v6 we have no positive detection (at the 3-ct confidence level) 
of the source and a 3-a upper limit is needed. We define the 3-cr upper limit as the number x that gives a probability to 
observe m or less counts equal to the formal 3-cr Gaussian probability, i.e.: 

P{<m,X + b)=PGauss{i<7) (1) 

Assuming Poisson statistics, eq. (1) becomes: 

^-i.+b) J2 I^±^ = 2.7 X 10-3 (2) 

In the limit of large numbers cq. (2) reduces to 

m = X + b — 3vx + b (3) 

By solving (3) with respect to x we have 

9 + 2to + 3V9 + 4to - 26 
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Table 1: Journal of BcppoSAX LEGS and MEGS Observations 



Date (UT) 




LEGS 






MEGS 






t"(s) 


Flux'' (x lo- 


-3 ctS S-1) 


^^(s) 


Flux'=(x lo- 


-3 CtS S-1) 






st 


S2 




st 


S2 


1998 Apr 26.334-27.458 


24483 


4.0 ±1.3" 


<4.0 


37220 


4.6 ±0.6 


2.4±0.5 


Apr 27.469-28.160 


13566 


<7.0 


<7.0 


21805 


4.5 ±0.7 


<2.5 


May 02.605-03.621 


1016.5 


- 


- 


31975 


3.0 ±0.5 


1.4'=±0.5 


Nov 10.754-12.004 


16961 


<6.0 


< 6.0 


53122 


1.8 ±0.4 


< 2.0 



" On source exposure time 

'' In the energy range 0.1-4 keV 

'^ In the energy range 1.6-10 keV 

"^ All uncertainties are at l-cr; upper limits are at 3-a 

'^ The 3-cr upper limit is 1.9 x 10^^ ^ts s~^ 



Table 2: Fits to the BeppoSAX LEGS+MEGS Average Spectrum of SI in April-May 1998 



Model 



ai 



a2 



^break 

keV 



kTTHB 

keV 



kTsB 
eV 



Red. 



X 



^ D.o.F. 



Single Power-Law 
Thermal Bremsstrahlung 

Broken Power-Law 

Power-Law + Black Body 

Th. Bremss. + Black Body 



1.0 ±0.3 


- 


- 


- 


- 


1.18 


12 


2.0 ±0.5 


0.7 ±0.4 


1.4 ±0.6 


7-^-1%' 


- 


1.66 
0.49 


12 
9 


0.5 ±0.3 


- 


- 


- 


90t^Ji 


0.49 


10 


- 


- 


- 


ifi + lSO 


90tfo 


0.46 


10 



° Uncertainties are at 90% (I.6-17) confidence level 



10 

Fig. A1.— BeppoSAX WFC (top) and GRBM (bottom) light curves of GRB980425. The onset of the GRB, indicated by the zero 
abscissa, corresponds to 1998 April 25.909097 (i.e., 5 seconds earlier than the GRBM trigger time). The vertical bars represent the typical 
1-(T uncertainty associated with the individual flux points. 

Fig. A2. — Discrete Correlation Function between WFC and GRBM light curves. The inaxirnuin at ~5 seconds indicates positive correlation 
between the two curves with the WFC light curve lagging the GRBM one by that temporal lag. 

Fig. A3.— BeppoSAX MFCS images of the field of GRB980425 in the energy band 1.6-10 keV referring to the epochs: (a) 1998 April 
26.334-27.458; (b) April 27.469-28.160; (c) May 2.605-3.621; (d) November 10.754-12.004. The local background has been subtracted and 
a smoothing has been applied with a window function of width comparable to the MFCS detector point spread function (3' at half power 
diameter). The source SI is clearly seen in all images. The fainter source 82 is detected in the first part of the first pointing, panel (a), and 
in the May 1998 pointing, panel (c). Note that in May 1998, panel (c), the pointing center was significantly displaced with respect to the 
April pointing, which accounts for the presence of the bright source toward North-East, about 9' from the image center, not visible in panels 
(a) and (b). In each image, the most external contours represent 2-a fiux levels. 

Fig. A4. — BeppoSAX LECS and MECS spectra of SI fitted with (a) a single power-law and (b) a power-law plus black body. The lower 
panels show the ratios between the data and the model. See Table 2 for the fit parameters. 

Fig. A5. — BeppoSAX LECS and MECS spectra of SI fitted with (a) a thermal bremsstrahlung and (b) a thermal bremsstrahlung plus 
black body. The lower panels show the ratios between the data and the model. See Table 2 for the fit parameters. 

Fig. A6. — BeppoSAX LECS and MECS spectra of SI fitted with a broken power-law. The lower panel shows the ratios between the data 
and the model. See Table 2 for the fit parameters. 

Fig. A7. — BeppoSAX MECS light curves in the 2-10 keV band of the X-ray sources SI (open squares) and S2 (filled circles) detected in 
the GRB980425 field. The WFC early measurements in the same band are also shown (stars). Uncertainties for the NFI measurements are 
1-(T. The \-a error bars for the WFC points are smaller than the symbol size and have not been reported, (b) Same as (a) for source SI 
only. The fits to the temporal decay with a power-law of index ~ 0.2 and with an exponential law of e-folding time -^ 500 days are shown 
as dotted and dashed curves, respectively, (c) Same as (a) for source S2 only. The first NFI measurement of S2 in (a) is replaced here by 5 
points obtained by integrating and averaging the flux in shorter time intervals. The dotted lines represent the power-laws of indices p ~ 1.6 
and p ~ 1.4 connecting the last WFC measurement and the first and last of the 5 NFI points of the April 26-27 light curve, respectively. The 
extrapolations of the power-laws to the time of the third NFI observation (May 1998) fall below the lower bound of the data point, although 
they are compatible with it at the ^ 2.5(t level. 
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Fig. 1 - BeppoSAX WFC (top) and GRBM (bottom) light curves of GRB980425. The onset of the GRB, indicated by 
the zero abscissa, corresponds to 1998 April 25.909097 (i.e., 5 seconds earlier than the GRBM trigger time). The vertical 
bars represent the typical 1-a uncertainty associated with the individual flux points. 
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Fig. 2 - Discrete Correlation Function between WFC and GRBM light curves. The maximum at ~5 seconds indicates 
positive correlation between the two curves with the WFC light curve lagging the GRBM one by that temporal lag. 
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Fig. 4 - BeppoSAX LEGS and MEGS spectra of SI fitted with (a) a single power-law and (b) a power-law plus black 
body. The lower panels show the ratios between the data and the model. See Table 2 for the fit parameters. 
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Fig. 5 - BeppoSAX LEGS and MEGS spectra of SI fitted with (a) a thermal bremsstrahlung and (b) a thermal 
bremsstrahlung plus black body. The lower panels show the ratios between the data and the model. See Table 2 for 
the fit parameters. 
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Fig. 6 - BeppoSAX LEGS and MEGS spectra of SI fitted with a broken power-law. The lower panel shows the ratios 
between the data and the model. See Table 2 for the fit parameters. 
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Fig. 7 - BeppoSAX MECS light curves in the 2-10 keV band of the X-ray sources SI (open squares) and S2 (fiUed circles) 
detected in the GRB980425 field. The WFC early measurements in the same band are also shown (stars). Uncertainties 
for the NFI measurements are 1-cr. The 1-cr error bars for the WFC points are smaller than the symbol size and have not 
been reported, (b) Same as (a) for source SI only. The fits to the temporal decay with a power-law of index ^ 0.2 and with 
an exponential law of e-folding time ^ 500 days are shown as dotted and dashed curves, respectively, (c) Same as (a) for 
source S2 only. The first NFI measurement of S2 in (a) is replaced here by 5 points obtained by integrating and averaging 
the flux in shorter time intervals. The dotted lines represent the power-laws of indices p ~ 1.6 and p ~ 1.4 connecting 
the last WFC measurement and the first and last of the 5 NFI points of the April 26-27 light curve, respectively. The 
extrapolations of the power-laws to the time of the third NFI observation (May 1998) fall below the lower bound of the 
data point, although they are compatible with it at the ^ 2.5a level. 
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